Detection of gravity waves by phase modulation of the hght from a distant star 



b 



in 
o 
o 

(N 



o 



G. B. Lesovik'", A. V. Lebedev"'^ V. Mounutcharyan", T. Martin'' 
" Landau Institute for Theoretical Physics RAS, 117940 Moscow, Russia and 
Centre de Physique Theorique et Universite de la Mediterranee, Case 907, 13288 Marseille, France 



We propose a novel method for detecting gravitational waves (GW), where a light signal emitted 
from a distant star interacts with a local (also distant) GW source and travels towards the Earth, 
where it is detected. While traveling in the field of the GW, the light acquires specific phase modu- 
lation (which we account in the eikonal approximation). This phase modulation can be considered 
as a coherent spreading of the given initial photons energy over a set of satellite lines, spaced at 
the frequency of GW (from quantum point of view it is multi-graviton absorption and emission 
processes). This coherent state of photons with the energy distributed among the set of equidistant 
Ch ' lines, can be analyzed and identified on Earth either by passing the signal through a Fabry-Perot 

^ , filter or by monitoring the intensity-intensity correlations at different times. 



^ . INTRODUCTION 



J> ■ The detection of gravity waves (GW) has stimulated a lot of interest for decades. There are two major GW 

', detection concepts: acoustic and interferometric detection. The acoustic method deals with a resonance response of 

' massive elastic bodies on GW excitations. Historically the acoustic method was proposed first by J. Weber 0] where 
he suggested to use long and narrow elastic cylinders as GW antennas. Although a significant progress has been 
achieved in fabrication and increasing sensitivity of such type of detectors ,3;, ^] the interpretation of obtained data 

\ is still far to claim undoubtedly the detection of GW. On the other hand a considerable attention has been shifted 

• recently to more promising interferometric detection methods. The interferometric gravitational- wave detector like 
' Laser Interferometric Gravitational Wave Observatory (LIGO) and VIRGO represents a Michelson interferometer 

^ I with a laser beam split between two perpendicular arms of interferometer. The principles of operation of such type 

Q \ of detectors are reviewed in Refs 0; M IS liol lilt . The action of gravitational waves on an interferometer can be 

^ ■ presented as relative deformation of both interferometer arms. A gravitational wave with dimensionless amplitude h 

\ induces the opposite length changes 51/1 — l/2hcosQ.t in each arm of the Michelson interferometer, where I is the 

• length of the arm, Q is the gravitational wave frequency. These length changes produce opposite phase shifts between 
^ ' two light beams in interferometer arms, when interference occurs at the beam splitter of Michelson interferometer. 

■ The resulting phase shift of a single beam of light spending time r in the interferometer can be written as JTj 

where lu is the light frequency. This phase shift results an intensity signal change of the light from interferometer 
beam splitter hitting the photodctcctor. 

The main problem of the acoustic and interferometric methods that they both deal with gravitational waves with 
extremely small amplitudes of the order h ~ 10"^^ [l^l reached the Earth from deep space. One can see from 
eq. that for gravity wave frequencies in the 1 kHz range, f2 ~ 10^ Hz, and for the light in visual frequency range, 
u! ^ 10^** Hz, one has the maximum phase shift of the order (50 ~ 10"^" for interferometer arms length of the order 
150 kilometers. Such extraordinarily weak effect requires a exceedingly high detector sensitivity both acoustic and 
interferometric detectors. 

Alternatively, GW detection may be based on effects associated with propagation of light or electromagnetic waves 
in gravitational fields. There are two primary effects for the light in constant gravitational field i) the deflection of 
light rays near massive bodies 'iS*! and ii) the Shapiro effect accounting for integrated time delay of the signal passing 
near a strong source of gravitational field [l^ . The same effects have to be observed for light propagating in gravity 
waves: the gravitational waves have to induce a weak time dependent deflection of light ray propagating through these 
waves and also lead to gravity-wave-induced variation in time delay. The idea to use astrometry to detect periodic 
variation in apparent angular separation of appropriate light sources was proposed by Fakir |l5j . It is shown that for 
a gravity wave source located between the Earth and the light source (with line of sight close aligned to gravity wave 
source) a periodical variation of the order Acj) ~ 7r/i(A) in the angular position of the light source has to be observed. 
Here h{A) is the dimensionless strength of the gravity wave at distances of the order gravitational wavelength A which 
is many orders of magnitude greater than the strength of the same waves when they reach the Earth. On the other 
hand one can directly measure the variation of the integrated time delay induced by gravity waves on the light emitted 
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FIG. 1: A star emits light which is deformed by the gravity field, but its phase is modulated, an effect which can be detected 
on the Earth. 



by distant star which passing through space region with strong gravity wave. The idea to use timing observation for 
detection of the gravitational waves was suggested first by Sazhin llol and then this problem was studied in details by 
several authors llsL [iQj . The estimations carried out in Ref. [17| give the following answer for the rate of change 
in the gravity-wave-induced time delay, f: 

\r\^\h{r^D)\ (2) 

where h{r = D) is the gravity wave strength at distances of the order of impact parameter D for the light beam 
passing near the gravity wave source. To put some numbers let us consider a dissymmetric rotating neutron star 
with spin frequency of the order 10'^ Hz and gravity wave amplitude of the order H ~ 10~^ cm then eq. 10) gives 
|f| ^ 10^^^ while the same effect measured near Earth results |f| ^ 10^^^ assuming that the neutron star to be at a 
typical distance of a few kiloparsecs from the Earth. One can see that all effects due to gravitational waves near the 
gravity wave source are several orders of magnitude stronger than the same ones on the Earth. 

In the present paper we want to suggest a new gravitational waves detection method based on the interaction of 
the photon with gravitational waves. Assuming that the photons from distant star passing near gravitational wave 
source, where the photon-gravity wave interaction assumed to be strong, the photon-gravity wave interaction leads 
to relatively strong modulation in time of photon frequency. The latter allows to analyze this modulation of the 
photons reaching the Earth by means of standard optical methods including the Fabry-Perot analysis and quantum 
photon correlations measurements. It is important that while the interaction of photons with gravitational wave 
is rather weak (proportional to strength of gravitational wave) the frequency modulation can be accumulated over 
large distances during the photon propagation that could result in an experimentally measurable effect on the Earth. 
In some aspects our treatment of the photon-gravity wave interaction resembles the effect of photon acceleration 
by gravitational wave p(3|, where the photon propagating in plane gravitational wave long enough time acquires a 
considerable increase of the frequency. However while in Ref. j2fl| the photon-gravity wave interaction was treated in 
the frame of reference of the photon, we have considered this interaction in the point of observation frame coordinate 
which seems to give rather the widening of the initial monochromatic photon wave packet than the increase of the 
photon frequency. 



PROPAGATION OF LIGHT NEAR THE LOCALIZED SOURCE OF GRAVITATIONAL WAVES. 



In this section we consider how the plane electromagnetic waves interact with gravitational wave field emitted by 
some localized source. The situation we have in mind is depicted in Fig. ^ The light signal originates from a distant 
star S and travels towards the Earth to be detected by an observer O, but along the way it interacts with a gravitational 
body M, which also emits gravity waves. We consider the trajectory of the light ray in the two dimensional plane 
formed by the star S, the body which emits the GW, Af, and the Earth. Provided that the wave length of the GW 
is large compared to the wave length of the light signal, we will use the eikonal or geometrical optics approximation 
to describe the interaction of the light in gravitational waves. The propagation of the electromagnetic waves in the 
field of gravitational waves was considered many times and wc can refer the reader to several papers dealing with this 
questions in the geometrical optics approximation _7,, .^L 22io 2^ 2Ae 2^ ■ Our subsequent analysis of this problem will 
be carried out in close analogy with the paper of Sazhin and Maslova [2^ , where they have considered the structure 
of electromagnetic field in Fabry-Perot resonator in the field of gravitational waves. However in the present paper 
we will be interested in a different aspect of the interaction of the light beam passing in immediate proximity to the 
gravitational wave source, where the amplitude of the gravitational waves is assumed to be strong in comparison with 
the amplitude of the same waves reaching the Earth. It is important that since the light interacts with gravitational 
waves in the region where they are strong, the resulting photons frequency modulation seems to be appreciable to 
detect it on the Earth. 



The eikonal equation. 



In the geometrical optics approximation the propagation of the hght ray is described by the eikonal equation: 

where g^^ ~ g^^^ — h^^ the metric tensor associated with the Schwarzschildian static metric g'^^ of the object M, which 
is perturbed by the tensor K^^ associated with the GW emitted from M • Here we will neglect the true form of the 
static metric g}}^ assuming that light ray propagates in the flat Minkovsky space g}^ = rf^ = diag{l, — 1,— 1,— 1}. It 
is known that the static metric leads both to the deflection of the light ray near the massive body and time delay of 
the light signal passing along the ray (so called Shapiro effect). It can be rigorously shown in the subsequent analysis 
that the former effect of the light deviation leads to the negligible correction to the additional phase accumulated by 
photons due to the interaction with GW and one can assume the light trajectory to be a straight line. As for static 
Shapiro time delay it can be easily incorporated to the final answer for photon phase and actually does not affect on 
the our subsequent detection method of modulated photons. 

Let "00 = —^^t + ky be the eikonal in the absence of gravitational waves perturbation with metric if^ (we assumed 
the light ray to be propagating along the y axis in the plane formed by three bodies) . Assuming that the perturbation 
is small, in the presence of gravitational waves the eikonal becomes V' = V'o + "^ij where V'l a small addition to the 
eikonal V'o computed to first order in - it satisfies the equation: 

dt + dy 2 ^' ^ ^ 

where we have introduced the notation F{t, r) = h^^ + K^^y — 2h^y . Using the Green's function formalism one can find 
the general solution of this equation. Then to first order in the gravitational wave perturbation, h^^ , the additional 
phase acquired by the photon due to the interaction with the GW, can be written for arbitrary time dependence of 
the perturbed metric: 

V 

1 ^ [ Tj^f^ y~y' 



Mt,y) = -^- I F{t-^^-^,y'] dy' (5) 



The emission of GW. 

Assuming that the GW tensor, ft,*'^, is a small perturbation to the static flat metric in the quadrupole approximation 
one can write the coordinate components of the GW tensor as : 

~ 2G Qc.p{t-r/c) 

= ^ (6) 

where G is the Newtonian gravitation constant, hi^ = — ^rjikh (where h = h^), Qap is the quadrupole moment: 

Qo,p{t) = / foo{t,r)x''x^ Sx, (7) 



where Tiu = Tik — ^rjikT, Tik is the energy-momentum tensor of the GW source, T = T^. Eq. (0 implies that the so 
called the harmonic gauge condition, dkh'H — has been chosen for the components of GW tensor hik- 

Using the gauge condition one can find the /iqo and ft-oa components of the GW tensor. To do this let us write the 
equation dkhf — separately for coordinate and time components: 

dhgp _ dhgQ dhpa _ dhpo 
Combining the eq. ® with the first eq. ||SJ) the h"'^ components of GW tensor can be written as: 
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while from this equation and the second equation JSJl one has: 



2G 



c2 dx'^dxf^ 



Qapit - r/c) 



(10) 



Using the expressions for the components of the GW tensor write the exphcit expression for F{t,r) = 

/jOO _|_ fj^yy _ 2}i^y in the r. h. s. of eq. Taking coordinate derivatives and keeping in mind that along the ray 
trajectory dxr — D/r, dyV = y/r, dzr = (where D is the impact parameter) the F(t,r) takes the form 



2G 
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2G 
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(11) 



where all quadrupole moments Qap are assumed to be the functions of the retarded variable t — r /c. 

Solution of the eikonal equation. 

For simplicity the source of GW is assumed to emit only one frequency 17. We consider the two following situations: 
either the GW source has been emitting since t = — oo, or the signal is bounded in time, with an exponential decay. 
In the former case one can write Qa^it) = qafs sin fit. Then substituting the expression for F{t,r), see eq. into 
eq. (O one could write the solution for tpi as an integral over the variable z = y' /D — -^/l + {y' /Dy-. 
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where the variable = t — y / c accounts for the effect of retardation of the light signal during its propagation to the 
point of observation, Hap = 2Gqap/c'^, A is the wavelength of gravitational waves. Assuming that at the point of 
observation y ^ D, one can safely set the lower limit for all integrals in this expression equal to zero. Then for large 
values of the impact parameter 2itD /A 3> 1 one can asymptotically estimate these integrals and the main contribution 
to ipi is given by the following expression: 



i-M) - ^ 



AV Lu 



sin 



(13) 
(14) 



where in the last line the answer is rewritten in terms of the dimensionless strength of gravitational waves: h^x = 
HxxA~^. One could see that the main contribution to the photon phase ^ comes from the hxx component of 
gravitational radiation confirming the well known fact |l8Lll9| that only gravitational waves propagating perpendicular 
to the photon wave vector (directed along y axis in our geometry) have a considerable effect on the electromagnetic 
radiation. 

Let us now consider the case where gravity waves are emitted starting at i = decaying exponentially with 
characteristic time tq. In this case let us choose Qapit) to be equal to: 



lap{t) = qaf}d{t) exp ( --^ ) sinQt, 



(15) 
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where 9{x) = for a; < and 1 for a; > 0. Then for large values of the impact parameter, 2ttD/A ^ 1, the solution 
for ipi is given by the following formula: 

One can see that in this case at the point of observation the modulated contribution to the eikonal also decaying 
exponentially with the same characteristic time tq. 



The pre-exponential factor. 

Let us now go beyond the geometrical optics or eikonal approximation and calculate the pre-exponential factor for 
the photon wave function. The photon wave function ip{t, y) with momentum k satisfies the wave equation: 



1 



9^1^-.]-^^ (17) 



where g = detgik- Up to the first order contribution in K^^ one could find g = — (1 + /i) with h = h\. Taking the 
derivative over a;* and keeping only the first order terms in K^^ one could arrive to expression: 



dx^dx^ dx'^ dx^ \ 2 

One could see that under the harmonic gauge condition dihl. = with h}. = h\ — ^rjl.h the second term in above 
equation vanishes and the photon wave equation reduces to 

^ dx'^dx^ ( ^ 

Let us write the photon wave function as tf{t, y) — c{t, y)e^'^ , where ip is the eikonal of the photon interacting with 
gravitational waves. Then assuming that c = 1 + a, where a is the small contribution of the first order in h^'^, one 
could rewrite the above equation as 

it / d'^a dib da , , d'^ih , Bijj dib \ , , 

^ (a^ + + ^^oJ?^ - (^ + ) - (^«) 

The last term in this equation vanishes due to the eikonal equation ||2Jl. Since the free eikonal ipo = ky — wt is the 
linear function of the all second derivatives of ipQ are vanish and the remaining wave equation up to the first order 
terms in K^^ takes the form: 

^^kJ^^<2^^^k^^^^^^^^k_^^^Q^ (21) 

dx^dx'^ dx^ dx^ dx^dx'^ 

Assuming that the a(t, y) = a(t — y/c) thus the first term above vanishes one could finally find 

da da _ I / I d'^ipi d^ipA 

dt dy~ 2 uj dt^ dy^ ) ' ^ ' 

Substituting here the general solution for -01 (see eq. (O) one can express the r. h. s. of eq. H22|l through the function 
F{t,y): 

da da _ c f dF 1 dF\ 

m ""d^ " I y&iJ ^ c~dt J ( ' 

According to our previous analysis the main contribution to a{t, y) is given by the H^x component of gravitational 
radiation: 
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Substituting this expression to the eq. the solution for a(^) for large values of impact parameter 2t:D/K ^ 1 has 
the form: 



[^j cosine (25) 

Combining together eq. (|14|l with eq. (|25|) one can finally write the photon wave function at the point of observation 
in the form: 

g-j:w(C+rgSin05) 

'^'-^^)= l + r,17cosl]r 
where we have introduced the time Tg associated with integrated Shapiro time delay (see below): 

The same solution eq. (|26|l for the photon wave function is valid for the case where gravitational radiation are bounded 
exponentially in time (see eq. (|15|l ) however in this case Tg depends on time and equals to 

-.(0 = -^ (^)\-'e{Oe-^'-'. (28) 



Modulation due to alternating Newtonian gravitational field. 

The observation and detection of the gravitational waves undoubtedly provides us a new astrophysical tool which 
can give us a deeper insight on the processes like the neutron star formation, the processes happening in close binaries 
at the final stage of it evolution, etc. . . . However the observation of gravitational waves gives us also an additional proof 
of the General Theory of Relativity itself and it seems reasonable to consider the effect of photon phase modulation 
in Newtonian theory of gravitation where the gravitational potential instantly adjusts to the current configuration of 
moving bodies. 

Consider for example an oscillating neutron star or a rotating double star system which in Newtonian theory of 
gravity causes an alternating Newtonian gravitational potential (j)(r,t). This alternating field can also lead to photon 
phase modulation and results in principle in the same effect like for the case of photon in the field of gravitational 
waves. The aim of this section is to compare the effect of photon phase modulation between these two theories. 

As the rotation or the oscillation frequency of the stars is much less than the frequency of the photon we can work 
in the same eikonal approximation as in the previous sections. The Newtonian gravitational field contributes only to 
the diagonal metric elements: = h" = 20/c^, where - is a gravitational potential at the point of observation. 
At large distances compared to the size of the star system we can leave only quadrupole term in the potential 

^{r,t)=GQo.0^. (29) 

Here r - is the distance from the star to the photon, the coordinate axis are chosen to be the same as before, so 
r — y^y^ + D^, m - the direction vectors of r: Hx — D/ y^y^ + D^ , Uy = y/\/D^ + y^, Uz = 0, Qap{t) is the 
alternating tensor of the mass quadrupole moment of the star. The values of Qafi are determined by the parameters 
of the star, like orientation of the angular velocity, eccentricity, etc. We will further discuss the general case, assuming 
only, that due to rotation or oscillation Qap{t) = ga/jsinSlt. 

The correction to the eikonal due to rotation is given by the same equation © as for the GW case solution of the 
eikonal equation, where in our case F{t, y) — h^^ + /i^^ = 2h°^ 

V 

*n = -^^ j h''[t^y^,y^ dy' (30) 

— oo 

Substituting the expression for (j) into the formula above and making the substitution z — y' / D one obtains: 

, y/D 

= 2c(,2^j J (TT^^ sm^l7e + 2.-zJd., (31) 
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where we have mtroduced Hap = 2Gqafj/c'^- 

Then for large values of the impact parameter 2ttD/A 1 and large distances y/D ^ 1 one can obtain the 
following approximate expression for 

" V^'^"'"''^^ " Hyy)^^^^^ + 2H^yC0snC]. (32) 

One can see that the photon phase modulation due to the Newtonian alternating potential vanishes exponentially 
for large impact parameters in contrast to the GW modulation which vanishes as (see eq. p4f) '). Technically 

the difference between Newtonian and GW cases follows from the different y' dependence of the F{t, y') function. For 
the properly retarded GW radiation the solution for ipi, see eq. (I14f) . involves the oscillating integral over sinri[^ + 
{y' — y/Jj/y^ + ly^) / c], where the oscillating dependence on y' near the source \y'\ ~ D almost vanishes making the 
light-GW interaction more effective. Oppositely for the instantly adjusted Newtonian radiation one has in eq. H3()|l 
fast oscillating integral over y' with sinri(^ + y'/c) near vicinity of the gravitational wave source, which leads to a 
very small value of modulation. 



The Shapiro effect in GW. 

Using the results of the previous sections let us calculate the integrated time delay of the light signal propagating 
in the field of gravitational waves. Consider an arbitrary wave packet emitted by some distant light source traveling 
toward to Earth near a source of GW: 

^(0^/y'c.e— «^ (33) 

where ^ = t — y/c. Using the results for monochromatic, photons propagating in the field of gravitational radiation 
one can find the signal at the point of observation: 

= ip{^ + Tg smn£_). (34) 

One can see that the interaction with the gravitational radiation leads to alternating time delay of the signal St ~ 
TgSmflt, (note that the Tg is negative, see eq. (|27|l '). This effect was suggested independently first by Sazhin and 
Detweiler for detection of gravitational waves using the timing observation of a pulsar the line of sight to which 
passes near the source of GW. 



DETECTION OF THE GW 



Let us now turn to the discussion about how the photons modulated by the interaction with the gravitational waves 
can be measured in some realistic setup. 

Consider first the simple case where the modulated photons (|26|l simply hit a photodetector which can react at all 
photon's frequencies. However, any real photodetector reacts not on the vector potential A{t) of the photons in the 
photodetector directly, but rather on the electrical field, E(t), induced by the photons in this photodetector. Then 
the one photon photo-detection probability, P(t)At, to observe a photon during time interval At can be expressed 
as H^: 

- " ^ (e^-^ it. r) • E(+) (t, r) ) (35) 

where S is the area of the photodetector, Edec some characteristic energy describing the interaction of the photode- 
tector with incoming photons, a is the quantum efficiency of the photo-detection, &'^'>(t, r) is the positive frequency 
part of the electrical field operator, E'^~)(i,r) is the Hermitian conjug ate of E(+)(i,r): 
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At the point of observation we set r = to shorten the notations. Taking the time derivative of the photon wave 
function H26|) one gets the photo-detection probabihty equals to 



Aira ScAt ■^-^ hujk 

EdetL^ 



P{t)At = ^ Zj^'^fc^ 



ks 

^ s . f , huj , . ,„„. 

= 2a^At duj - — n(uj) , (37) 
r'' J Edet 

where r is the distance between the hght source and the Earth, n{Lo) duj is the number of the photons with frequency 
CO, emitted by the distant hght source per unit time, into a unit angular domain. 

One can see that the photo-detection probability does not feel any photon phase modulation due to GW and this 
probability is the same as it would be for the monochromatic photons emitted by a distant light source. Thus for such 
a simple setup it is not possible to retrieve any information about the gravity wave from the photodetector signal. 

Consider however a more realistic setup where the photodetector has a sensitivity edge lUs which means that 
photodetector does not react on the photons with frequency less than LUg. One can treat this situation as if the 
photons passed through a filter with transparency t(Lu) — 9[ijj~uJs), before hitting the photodetector. It is convenient 
to write this transparency as an operator acting on the photon wave function: t — 9[idt—LiJs\ ■ Then after passing such a 
filter the photon has a wave function equals to tipk{t) ~ 6[LUk{t) ~ LUs]ipk{t), where Ukit) — ujk[l + cos fit] — ujkf{t). 
Then for the photo-detection probability one has: 

P(t)At « 2a / dtj— 9\iuf(t) - LUs] n(uj). (38) 

J Edet 

Consider now the case where photons coming to the photodetector from a sharp Lorentzian spectral line with maxi- 
mum, 0^0, near the photodetector sensitivity edge luq « and width Tq: 

then according to the eq. (|38|l the photo-detection probability for the modulated photons H26|l has an appreciable 
contribution which is periodic in time at the GW frequency, and which is proportional to the strength of gravitational 
waves: 

P{t)Atfii PoAt[l + nTg— cosnt] , (40) 

where P^At is the photo-detection probability for non modulated photons. 

From the above analysis one can conclude that in order to detect the modulation of the photons by GW one should 
place before the photodetector some filter which has a finite frequency band. Then the formula for photo-detection 
probability (|38|) with a sensitivity edge has a rather general sense. In fact let us consider an arbitrary filter placed 
before the photodetector with transparency T{io). Then provided that the time spent by the photon in this filter (that 
is an inverse frequency band of the filter) is much smaller than frequency of the GW, 17, one can treat the modulated 
photons as monochromatic photons with slowly varying frequency Wfe(i) = lo-^SX ~\- rjrgcosfii] — LOkfit). In this case 
the photo-detection probability can be written simply as 

P{t)At » 2a ^ / dto^ T{Luf{t)) n{Lj) . (41) 
J Edet 

From this expression one can see that there are two diS'erent regimes of GW detection in this proposed setup. 
Consider the first regime where the distribution function of the light source, n{uj), has a sharp peak of the width Fq 
(or in more general case some sufficient irregularity like the edge of the spectrum) centered near the transmission 
window of the filter. Then, provided that the amplitude of the photon's frequency modulation, ujk^Tg, is much bigger 
than the width of the spectrum irregularity, uj^^Tg ^ Tq, the resulting photo-detection signal will be a periodic 
sequence of peaks. 

In the opposite case where the light spectrum n{uj) is a slowly varying function of the frequency within the trans- 
parency window of the filter one can consider n^uj) as a constant in eq. H41II . Then the time dependence of the 
photo-detection probability can be written as: 

P{t)At « 2a ^ ^ I dc.^^ Tic.) niu;) . (42) 
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This result can be understood as the time modulation of the number of incoming photons per unit time within a finite 
frequency interval. 

From another point of view one can note that the modulated photon wave function (|26() can be written as a 
superposition of waves with energies shifted by nfl (n integer): 

1 +00 

^-(^)- l + ^7.,cosllt E^"^-^)-"'^""""" 

^ —00 

where Jn is the Bessel function. In order to detect a gravity wave signal it is therefore necessary to provoke an 
interference between these Fourier components. To do it one can both investigate the modulated light signal passing 
through the Fabry-Perot interferometer or from the other hand study the intensity-intensity correlation of the photons 
coming on the photodetector. 



Analysis of the light signal with an interferometer 

In this section we consider in detail the setup where at the point of observation - the Earth - the modulated 
light signal before hitting the photodetector passes through an interferometer, which is characterized by a complex 
transmission amplitude t{uj). We will assume that the light passing through a Fabry-Perot filter comes from a single 
Lorentzian spectral line of the type (|39(l . For simplicity, let us first assume that there is only one FP resonance within 
this spectral line. The transparency of FP is written in the usual way: 

t{u;) = '^/^ , . (44) 

Here a derivation of the photo-detection probability is provided. It relies on the assumption U,Tg <^ 1, which is relevant 
for experimental situations (see below). 

Consider the propagation of a photon wave-packet ipk{t) (|26l) through the FP interferometer. The resulting wave 
packet after the FP can be written (i) = tip^. (t) , where t is the transparency operator of the FP in real time 
representation: 

i=Y[i9t-LOo+ir/2]-\ (45) 

where F is the width of FP resonance. However the real FP filter does not act directly on the vector potential of the 
incoming photons but rather on the electromagnetic field induced by the these photons. Using similar arguments the 
positive frequency part of electrical field operator E(+^(t, r), after passing the Fabry-Perot filter, can be written in 
terms of the transmission operator of the filter: 

E(+) it,r)^-^Yl \/5 ^)] ""^^ ■ (46) 

With these notations the photo-detection probabihty can be written as 

P{i) = 2a— I [iifkit)] * [i'Pkit)] n{ij) duj. (47) 

The action of the operator t on the fmiction i^kit) can be found using the Green's function formalism: 

t 

i^k{t)^\ J e-^-°(*-^)-*(*-^Vfe(r)dT. (48) 

—00 

Substituting this expression into eq. 147|l in the limit loq ^ T one finally obtains the following expression for the 
photo-detection probability: 

2 * * 

P(t) = 2a4^^ [n{cu)e-^'clw [ [ e-'^°^^-'^+i^^'+'^^l{T)^k{s) drds (49) 
r-^ Edet 47 J J 
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Performing the integral over uj and introducing the relative time £^ — t ~ s, and the total time t] — {t + s)/2, in the 
working limit flTg <^ 1 one can write the photo-detection probability as: 

t 2{t-ri) 

r-^ Edet 47 J 

Performing the integration over ^ in the limit Fg ^ F the leading contribution to the photo-detection probability is 
given by the following expression: 

00 


where we have made the substitution z = {t ~ rj) and we introduced the notation ro(t) = Fo|l + Qrg cosilij. 



"Strong" GW or narrow spectrum. 

Let us consider the limit where the frequency broadening of the initially monochromatic photon wave packet due 
to the interaction with gravitational waves, ST = LUoQTg, is much bigger than the width of the spectral line, Fq. The 
condition ST ^ Fq can be achieved both in the case where the amplitude of GW is strong enough or when one has a 
very narrow spectral line. In this limit one can safely neglect the dependence on |1 -I- flTg cosilij and put Fo(i) = Fq 
in eq. H51|l assuming flTg <C 1. Then the leading contribution to the photo-detection probability can be written in 
terms of the occupation number, n{uj), of the spectral line: 

Pit) = 4 — r' / n[ujo{t - z)] e-^'- dz (52) 

Edet J 


where u>o(t) = ujq(1 -\- flTg cos fit). This result has a very simple qualitative explanation which is consistent with our 
previous simple arguments for the photodetector with sensitivity edge (see eq. (|41f) '). Consider a photon wave packet 
coming to FP filter: ipk = C(i)e~*"'=(*+'^9*''"^*). Such a wave packet in the limit 17 <C Wfc can be interpreted as a 
photon with a frequency which is slowly varying in time ujk{t) = Wfe ( 1 -|- ^iTg cos nt). If however the time spent by this 
photon in the FP filter, F^^, much smaller than the characteristic time of the frequency change fZ"^ one can assume 
the photon at any instant of time to be monochromatic with frequency u!k{t). Then one can take advantage of our 
eq. (|41|l . Performing the integration over ui in eq. H41|l with Lorentzian spectral line n{Lu) 13911 one has 

P{t)^aTi-—^n[iVo{t)]T. 53 
r-^ Edet 

On the other hand this result can be derived directly from our elaborate analysis by putting n(a;o(t — z)) = n{uja[t)) 
in eq. (|52|l . The exponential factor in eq. (|52|l describes a time delay of the photon wave packet propagating through 
the FP filter. 

The result eg . 1153(1 describes the periodic sequence of symmetric peaks with a half period of gravitational waves it /Q, 
and widths equal to 

r^^^^-\ (54) 

If however the width of these peaks becomes smaller than the time needed for the photon to penetrate through the 
FP filter: r ^ F^^ then the form of peaks becomes asymmetric and the more accurate expression (|52|l has to be used. 

It should be noted that the expression (|52|l is valid not only for a single spectral line (pTi^ but also for arbitrary 
distribution function for the light source, where there is a sufficient irregularity in n{u>) like, for example, at the edge 
of the spectrum. Indeed if one consider the arbitrary distribution function n{uj) which changes much at frequency 
scale Aw around then it follows from eq. H52|l that one again has a periodical sequence of peaks for photo-detection 
signal in regime 5T 3> Aw. In this sense the described effects (|52ll . (|53|l have nothing to do with the interference 
of the different components of the photon wave packet (|43|l and deal rather with the effect of the photon frequency 
modulation u}k{t) = Wfe(l -|- Org cos J7t) due to the interaction with the GW. If however the amplitude of the frequency 
modulation ^F becomes smaller than the characteristic scale of the spectral function n(w) (the case of a broad 
spectrum) the effects associated with frequency modulation completely disappear. 
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"Weak" GW or broad spectrum. 

Consider now the opposite limit to the previous case when the ampUtude of GW is smaU or when the width of the 
spectral line is large enough so that the amplitude of the frequency modulation, ST, is much smaller than the width 
of the spectral line, Tq. It should be noted that this situation is the most common situation in practice since the 
amplitude of GW is extremely small even in the vicinity of GW source. In the limit ST <C Tq one can safely expand 
the integral expression in cq. (|51|1 up to the lowest order in the small ratio ST/Tq ^ 1 and perform the integration 
over z. Keeping all terms of the type 1 + flTg cos ilt in the working assumption F ^ the photo-detection probability 
can be now written as: 

S Too ( \ 

P{t) « QfTT — n(wo)r [\-^Tg COS VLt - 2—2" cos 2m ) (55) 
r Edet \ / 

One can see that in the case of a broad spectrum one obtains small oscillations against a huge constant background 
with frequencies fl and 2fl suppressed by small factors flTg <C 1 and JF^/Fq ^ 1, correspondingly. However as one can 
see from eq. (|55|) even for the extreme case of infinitely broad spectral line still there is a time dependent contribution 
to the P{t): 

P(t) w ckTT 4 n{ujo)T [l - cos nt] (56) 

It should be noted that this expression does not depend absolutely on the form of the spectra of the distant light 
source and the time modulation of the photo-detection signal defined only by the strength of the gravitational waves. 
The amplitude of the photo-detection probability in this case proportional to the number of photons n(wo)F coming 
to the photodetector per unit time and practically one can use the FP filter with appropriate width F to collect the 
enough number of photons. 

The second important point which has to be emphasized here is that the time dependent term in eq. (|56|l has 
nothing to do with the effect of the frequency modulation mentioned above (see eq. 153|) and presents an absolutely 
different effect associated with interference of the different components of the photon's wave packet (|43|l after passing 
the FP filter. To clarify the nature of this interference effect we will consider in the next section an example where 
there are two extremely narrow FP resonances, F ^ f2, within the broad spectrum of the distant light source. 



Two resonances case. 



Let us consider a more general case where two narrow resonances at frequencies ± A/2 are present within the 
spectral line (A is the separation between resonances). In this case the transmission amplitude of the FP filter, t{uj) 
again can be written as an operator acting on the e.m.f. induced by modulated photons: t = ti + t2, where ti(2) 
corresponds to the first (second) resonance: 



1(2) 



idt 



A T 



(57) 



where 2(j) \s the relative phase difference between resonances. 

Using similar type of arguments and the same method of calculation as for the case of FP filter with a single 
resonance, in the limit Fq S> fi, Fq S> F one arrives at the following expression for the photo-detection probability: 



Pit) =a^^n{uj^)To^ [Pi (0 + 2Pi2 (t) +P2 {t)] 



(58) 



where 



Pm{t) = 

Pl2{t) = 



To{t-z)e 



-Tz 



dz 



(JFcosrj(t-z)± A)2 + £l(izf) 

To{t - z) cos(Az - (j)) e"^^ dz 
ST^ cos2 n{t -z) + 



(59) 
(60) 
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Consider now the extreme limit of infinitely broad spectrum: Fq — » oo, then in the limit F ^ F ^ A one has: 

P{t) = air— n(wo) (2F) [l - Qrg cos nt] . (61) 

Comparing this result with the expression for the photo-detection probability for the single resonance case (|56|l one 
can see that in the regime of a broad spectrum Fq — > cxd the photo-detection probability does not depend on the 
number of resonances in FP filter and depends rather on the total transparency of the filter (2F in present case). 

Consider however the opposite limit of extremely narrow resonances F <C ^2 (this is hardly possible in practice) . 
Then the resonance factor with respect to separation between resonances, A, appears in the expression for photo- 
detection probability: 



, , S TiLOn 

P(<)«a^-— ^n(cc>o)(2F) 



^ f^TgF Fcos(rJi-(^) (^7-A)sin(m- 



F2 + (fi-A)2 



(62) 



One can see that there is a parametric resonance for the photo-detection signal at A = il. The photo-detection 
probability at the resonance equals to: 

/'a=j2(0~« — Tr^"(^o)(2F)[l-— ^ cos(r!t-0)], (63) 

while for |A — fij ^ F (out of resonance) the term which depends on time in the photo-detection probability is 
saturated by the small factor F/|A — ilj. 

A similar resonance appears in regime F ^ for the case where the width of the spectrum n(ixi) is broad but 
(ci;o/Fo)^^lTg ^ 1 so one can neglect all factors |1 + firg cos51t| in eq. H58|) : 



' ^det 



5T^ Tcos{2nt-(j)) + {2n-/^)sm{2VLt- 



^ Fg ^ ■ F2 + (217 - A)2 



(64) 



One can see that the parametric resonance at A = 212 appears in this case. Now we argue that the resonances at A = 11 
and A = 217 H62I) , H64I) arises because of the interference of different components of the photon's wave packet H26|) . To 
show it let us simply write the transparency of the FP filter as t(Lo) — 5t{uj — loq) -\- (5(a; — — A) where 5y{lo) is a 
narrow peaked function with finite width F <C 17. Then using the expansion of the photon wave packet H43|) one can 
write the photo-detection probability in terms of the sums over different components of the photon's wave packet: 

P{t) = 2a— J2 / TT^ "(^) '^^ e*'"-'")"*i*(w + nn)t{uj + ml7) J„(tjrg) J™(wrg) (65) 

Let us assume that the distance between resonances are exactly equal to Nil — A where N is an integer. Then 
substituting the transmission amplitude t{uj)mto this expression and making the integration over lo one gets; 

P{t) w 4a4 f2 "('^o ~ ?t(^o - nfl) y2^^^^_^^y^^ ^ J„[(coo-nl7)rg] J„+jv[(t^o-ri^^)rg] cos Nilt] (66) 

^ ^det 
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Then for the special choice N = 2 assuming that the argument of the Bessel function is typically large one can put 
Jn+2{x) ~ ~Jn{x) for a; 3> 1. Then the photo-detection probability can be written as 

P(i) w4a^ (l-cos217i) ^ n(^o - nl7)F ^^""^^ """^ J^[(c^o - nl7)Tg]. (67) 



n— — oo 



One can see that the term which is alternating in time in P{t) arises due to interference of the different components 
of the photon wave packet (|26|l . 



Analysis of intensity correlations 



In the previous section, we showed that by appropriate filtering of the light signal, it is possible to extract oscillations 
associated with the past interaction with the GW. While this previous proposal is promising and attractive, one can 
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take an alternative route for the detection of the GW signah the measurement of intensity-intensity correlations of 
the photo-detection signal. Consider the correlator ((/(ti)/(t2))): because of the time modulation due to the presence 
of the GW, it no longer depends on the variable ti — t2 only. Using the definition the intensity operator, one can see 
qualitatively that the Fourier components of Eq. I|43|l will lead to a two particle interference effect. 

The quantity to be measured on the Earth is thus the intensity correlator of the photo-detection signal which is 
the time ordered product of electrical field operators: 

((: /(<i)i(i2) :» dhdh = «2^!£yMi2 ((E(-)(ti)E(-)(<2)E(+)(t2)E(+)(ii)» (68) 

Substituting here the expression for the electrical field operators H36(l and performing the quantum average one gets 

((: Hti)Ht2) ■.))dtidt2 = 2a2^^^^ / d^^du;2e-^"^[*i-*^+^«(^'""*i-^'""*=)l 



n{oji)n{u2) e»-i[*i-t2+r,(sinnti-si„m2)]^ 

El. 



(69) 



28^ dtidt2 I ^r, , , , . r^, . r-v, M|2 



This correlator can be rewritten as 

((: i(ti)i{t2) ■.))dtidt2 - 2a^ -^ - t2 + Tg(sinmi - sinm2)]r, (70) 

where 

/W = / ^n{u:)e-^-Uu; . (71) 

The function f(t) typically decays exponentially on a time scale given by the inverse of the line width Fq 3> fi. 
Thus one can safely linearize the dependence of the intensity-intensity correlator with respect to the relative time 
T = ti- t2: 

((:/(f + T)/(t):)) =2a2^|/[r(l + f]TgCOsm)]|' (72) 

where t — ti — t2 is the relative time and t is the total time. Let us now calculate the spectral power, S{t) of the 
photons' noise defined as 

S{t) = j{{:i{t + T)i{t):))dT (73) 

Substituting here the expression for photon intensity correlator one can show that 

S{t) = — —r / n (w) duj (74) 

One can see from this equation that the intensity correlation contain periodic oscillations - although with a small 
amplitude - in the total time. Note that contrary to the Fabry-Perot diagnosis, here a rather accurate measurement 
is implied, as the contribution to other noise sources (due to scattering...) needs to be minimized compared to this 
periodic signal. In particular, it requires a large time acquisition window to filter out spurious fluctuations. 



Conclusion 



Let us now summarize all results obtained through the article and give some estimations and limitations for the 
proposed gravitational waves detection method. The primary effect we have dealt with is the modulation of the 
photon phase due to interaction with gravitational field have been found in the geometrical optics approximation and 
described by eq. for the photon wave function: 



(75) 
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where ^ t — r j c and C(<^) = [l + f^Tg cosri^] ^ is the pre-exponential factor. This effect of phase modulation leads 
immediately to well known Shapiro effect of the time delay of the signal propagating near the source of gravitational 
field. For the case of the source of gravitational waves this time delay is an alternating function of time Tar = 
r/c — Tg sin Or, where Tar is an arrival time of the signal, r/c the traveling time of the signal in flat Minkovski space 
and 

'^-(^^ (^)\-\ (76) 

Here h^xi^) is the dimensionless strength of gravitational wave near the GW source, D impact parameter, — 2ttc/ K 
is the frequency of gravitational waves. The Shapiro effect in fact seems to be a very promising candidate for gravity 
wave detection experiment since the time delay (or the additional alternating in time phase) mainly accumulated near 
the source of gravitational waves where the strength of GW many order of magnitude greater when they reach the 
Earth. However in the proposed experiment E3 | it requires a pulsar on the line of sight which passes near the source 
of gravitational radiation. In this paper we propose a gravitational waves detection method based on the Fabry- 
Perot interference analysis (or equivalently time correlation measurement) of the light signal from arbitrarily light 
source (not necessarily a pulsar) passing near the strong source of gravitational radiation. Then the photodetector 
signal contains an alternating in time component (with GW frequency) proportional to the strength of gravitational 
radiation near the the GW source (see eq. EEl lfil|l : 

P{t) = Pa{l-VLTgCOsVLt), (77) 

with Po a- photodetector signal corresponding to a light coming to a FP filter from distant star. To estimate an 
effectiveness of the proposed method let us first estimate a brightness of the light source needed to resolve an alternating 
component in the photo-detection signal. Collecting the signal over long time Tobs allow to better resolve the alternating 
component. Assuming that during the observation time A^obs photons coming from a distant star hit the photodetector, 
in the limit of Poissonian statistics, {SN^]^^) — A'obs one has the following limiting requirement: 



r!rg(Z?) (78) 
in the extreme case where the light passes in immediate proximity {D w A) to the GW source one has the condition: 



/i(A)0Vo;^»l- (79) 

Then assuming that the distant star emits in the frequency band Alo and that the diameter of the telescope equals 
to d one has: 

where L is the brightness of the distant star, R is the distance to the Earth. 

The most promising candidates of gravitational waves sources appropriate for proposed detection method are the 
periodic sources of gravitational radiation such as close binaries at the final stage of evolution or a asymmetric rotating 
neutron stars, allowing to collect the signal long enough time to resolve the tiny gravitational wave modulation on 
the stochastic Poisson noise. Let us give some estimations, assuming the total mass, m of binary system of the order 
of the mass of the Sun, and rotating with a frequency v ~ 10^^ Hz. Then the dimensionless strength of gravitational 
waves at the distance r from the emitter, according to quadrupole approximation, is given by 

/z«4^im5/3:.2/3 (81) 
r 

assuming that the role of the distant star plays one of the components of the binary system r — [Gm/v^Y/^ one has 
h{r) A{Gmv)'^/^ / ~ 10~^. In order to resolve this modulation one needs about Aobs ^ 10^^ photons coming to 
a photodetector during Tobs- Assuming the the typical distance from the binary system to the Earth of the order 
1 kpa with d « 10 m, ~ 10^^ Hz and F/Acj ~ 1, Tobs ^ IQ^v'^ ~ lO^s one has the required brightness of the 
binary component, according to eq. H80|l of the order L ^ lO^is where Ls is the brightness of the Sun. Of course 
the requirement to have such a bright component of the binary system makes the observation of the proposed effect 
rather problematic in practice. 
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In reality the situation with binary systems is even worse due to the presence of the Doppler effect for the light 

emitted from one of the rotating components. In fact the Doppler effect results in the same frequency modulation of 
the emitted light as the gravitational waves modulation and the amplitude of the frequency modulation due to the 
Doppler effect is of the order of 



For parameters described above it results in an effect of the order 10^ that is five orders of magnitude greater than 
the frequency modulation due to gravitational waves. 

The situation however could be much better for very close neutron star binaries at the last stage of evolution. In 
this case the rotational frequency can be of the order of 10^ Hz and results in a dimensionless strength of emitted 
gravitational waves of the order h{A) ^ 10~^. Assuming that the distance to the binary neutron star system is of 
the order 1 kpa and that there is a bright distant star situated on the line of sight to the binary from the Earth with 
impact parameter D ~ A and collecting the signal during the time Tobs ~ 10^v~^ ~ 10 s one has the following lower 
limit for the brightness of distant star L ~ lO^Lg- 
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